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We investigate the superconducting states that arise as a consequence of mesoscopic confinement 
and a multi-component order parameter in the Ginzburg-Landau model for p-wave superconductiv¬ 
ity. Conventional vortices, but also half-quantum vortices and skyrmions are found as the applied 
magnetic field and the anisotropy parameters of the Fermi surface are varied. The solutions are 
well differentiated by a topological charge that for skyrmions is given by the Hopf invariant and 
for vortices by the circulation of the superconducting velocity. We revealed several unique states 
combining vortices and skyrmions, their possible reconhguration with varied magnetic held, as well 
as the novel temporal and held-induced transitions between vortical and skyrmionic states. 

PACS numbers: 74.78.Na, 74.25.Ha, 74.25.Dw, 74.20.De 


I. INTRODUCTION 

Strontium ruthenate, Sr2Ru04, is according to the¬ 
oretical predictions the best candidate to date to host 
p -wave superconductivity. Generally speaking, the order 
parameter in superconductors describes the spatial pro¬ 
file of the gap function, A^/fc). The order parameter in 
p -wave superconductivity is an odd function of the wave 
vector k, unlike the s-wave superconductors where it is 
an even function of Following the notation of Balian 
and Werthamer, the p-wave order parameter read^^ 
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or in a short notation A(fc) = i (d(fc) • d) ay, where d{k) 
transforms as a vector under rotations and ai are Pauli 
matrices. Microscopic calculation of the superconduct¬ 
ing gap is a highly demanding task that requires a de¬ 
tailed knowledge of the pairing mechanism which in many 
cases is not available. What remains then is to exploit 
all the symmetries (continuous and discrete) exhibited by 
the material under consideration and build a model that 
will depend on certain number of parameters. The pos¬ 
sible superconducting order parameters that have been 
reported for p-wave superconductors required a detailed 
description of the crystal structure of the considered 
material.l^ In that respect, strontium ruthenate (SRO) 
is a layered perovskite with a crystal structure similar to 
the well known high-T^ superconductor (La,Sr)2Cu04, 
where oxygen ions at the corner s o f an octahedron sur¬ 
round the body-centered Ru ion.hJ^ The planar layers of 
Ru 02 are separated by Sr layers that stack along the 
highly symmetric axis c. The Fermi surface of strontium 
ruthenate contains three sheets arising from the binding 
of the Ru and O ions within the same layer. Bindings 
between the Ru02 layers are weak due to the long sep¬ 
aration of the interplanar RuOe octahedra. The Fermi 


sheets a and (3 are both one dimensional (ID), while the 7 
sheet is two dimensional ( 2 D). For sufficiently high fields 
along the c direction, the Zeeman splitting between elec¬ 
trons with opposite spins demands vector dik) to lie in 
the Ru 02 layer or the basal plane,^i.e dzik) = 0 . Thus, 
just two components of vector d(fc) are left unknown. A 
rigorous analysis beyond the scope of this work found 
that among the possible representations for d^ ± idy, 
there is one which is 2 D and consequently breaking time- 
reversal symmetry.12121 This representation is suitable for 
strontium ruthenate since muon-spin relaxation (/rSR) 
measurements detected spontaneous magnetization be- 
lo w i ts superconducting critical temperature, Tc = 1.5 
kEEI fitting well to what was predicted in a formalism 
where the superconducting order parameter has at least 
two components. This 2 D representation introduces an 
internal degree of freedom, the chirality, that has been 
widely discussed and reported in recent works.^^^ For 
example, the work of Chung et al. provided the param¬ 
eters needed to stabilize a vortex which is not coupled 
to the electromagnetic potential and which carries half 
of the vorticity of a conventional vortex -in other words, 
a half quantum vortex (HQV). Skyrmion states carry ing 
topological charge defined by the Hopf invariant,g^j-g 
another hallmark of chiral superconductivity, and were 
obtained as stable solutions in bulk p-wave supercon¬ 
ductors within the GL model.l^ 

In this work we report various full vortex (FV), 
HQV and skyrmion states as the fundamental solu¬ 
tions of the GL mo del for mesoscopic chiral p-wave 
superconductivity.!^ It is well known in conventional s- 
wave superconductivity that confinement can stabilize 
superconducting configurations which in bulk systems are 
energetically unfavorable or even unattainable, e.g. non- 
Abrikosov vortex lattices, or vortices with phase wind¬ 
ing (j) = 27 rri, with n > 1 (giant vortices).Similarly, 
in mesoscopic chiral p-wave superconductors HQVs have 
been predicted to exist, owing the reduction of their oth¬ 
erwise divergent energy to the low dimensionality of the 
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systemP We found these HQV states in multiple forms, 
but also FV and skyrmionic states and transitions be¬ 
tween them as a function of the external magnetic field 
applied perpendicularly to the sample. We employed the 
time-dependent theoretical formalism, which allowed us 
to observe novel temporal transitions as well, related to 
peculiar entry and arrangement of HQ Vs and their tem¬ 
poral transformations into other topologies. 

The paper is organized as follows. Sec. [IT] presents the 
theoretical formalism and our analytical analysis of the 
first GL equation and the superconducting current. The 
boundary conditions imposed on our equations are de¬ 
rived from the latter expression. Sec. |III| then summons 
our findings for the superconducting configurations com¬ 
posed of HQV, FV and skyrmion states, obtained at weak 
coupling and considering a cylindrical Fermi surface. The 
transitions between states of interest as a function of the 
magnetic field are discussed in Sec. |IV[ while the tem¬ 
poral transformations are shown in Sec. |Vj The effect 
of anisotropy on the topological, vortical and skyrmionic 
entities is analyzed in Sec. |VI[ Ou r findings ans conclu¬ 
sions are summarized in Sec. ■Mil 


II. THEORETICAL FORMALISM 


After the above brief general description of strontium 
ruthenate, in what follows we show the Ginzburg-Landau 
(GL) equations that the order parameter, ^ 
must satisfy. The order parameter has two components 
(is chiral) as a consequence of the 2 dimensional represen¬ 
tation (F^) of the tetragonal group Z 34 ?i.l^The expansion 
of the GL free energy density up to fourth order in 'ipx,y, 
that fulfills the group symmetries, reads 


^ = K {\Dx'ipx\‘^ + iL’yV'yP) + + iDytpxl'^) 

+ 2Re^k2Dxi^x{Dy'tfjy)* + ksDxi’yiDy'ipx)*'^ -a 

+ + P2{'tl’xi’y - '>Px'^*yf + /JsIV'xHV’yP , (2) 

where a, ki and jSi, with i = 1,2,3, are parameters that 
depend on the details of the Fermi surface of the mate¬ 
rial under consideration. K — Dx,y denote 

the components of the covariant derivative. The time- 
dependent Ginzburg-Landau (TDGL) equations, used in 
our numerical approach,^ are the set of coupled differen¬ 
tial equations for the superconducting order parameter, 
and the vector potential 
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where (p is the scalar electric potential, B is the magnetic 
induction, is the effective mass, D is the phenomeno¬ 
logical diffusion coefficient, and a the electrical conduc¬ 
tivity. For convenience we set h = 1 and ms = 1/2. 
The second GL equation [Eq. @] is discarded in this 
work since the diamagnetic effects of superconductors are 
vanishingly small for a thin (effectively 2D) mesoscopic 
geometry. We use the symmetric gauge for the vector 
potential, A = {f x H)/2, with the magnetic held (H) 
directed along z. The scalar electric potential is set to 
zero since neither charges nor external currents are con¬ 
sidered in this work. In dimensionless units, where dis¬ 
tance is scaled to the coherence length, ^ time 

c2 

to to = 45 j magnetic held to the upper bulk critical held 
Hc 2 = and the superconducting order parameter 


to A+ = 



the hrst TDGL equation becomes 
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where n± = '^{Dx i iDy), tp± = ipx A: iijjy and r = 
/ 32 // 3 i. a straightforward calculation reveals the follow¬ 
ing important result, [Dx,Dy] = iH, which leads the 
operators n± to satisfy the commutator: [n 4 _,n_] = 
H. The external magnetic held, being constant, can 
be factored out from the above commutators, leading 
to [n+,n_] = 1 , which dehnes the algebra behind the 
Landau levels; L[± = II±/'/H. This algebra is dehned 
through the following commutators: [ 7 V,n+] = —n+, 
[V,n_] = n_; where N = H+H- is the particle number 


operator. Within the weak-coupling limit and consider¬ 
ing a cylindrical Fermi surface (7 sheet), all the ki pa¬ 
rameters are equal to {vlvy)/{v‘^) = 1/3, where brackets 
() denote averaging over the Fermi surface.l^ For this 
case the hrst GL equation reads 


= x 

3 L 


D'^+Ula++Ula- 
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where a± = ± iay)l2, are pseudospin or chiral op¬ 

erators acting on the space span by 'ip±. Ignoring the 
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nonlinear terms (linearized case), it is straightforward to 
show that the superconducting order parameter must be 
of the form: ^ = {(j)M, 4 ’n- 2 )'^, where (j)N is the state 
corresponding to the Landau level Within the 

superconducting formalism the number N turns out to 
be the vorticity of the order parameter. Then, one con¬ 
cludes that for chiral p -wave superconductors there is a 
vorticity difference two between the components of the 
superconducting order parameter. The full GL equa¬ 
tions, i.e. the linearized equation plus the nonlinear 
terms, are a complicated set of partial differential equa¬ 


tions with unknown analytical solutions. Therefore, in 
this work we solve this problem numerically. Due to the 
mesoscopic dimension of the sample under consideration, 
proper boundary conditions must be incorporated in the 
GL equations in order to pose the problem well. In what 
follows, the superconducting current is calculated for the 
general case, which includes the specific case where all 
kiS are equal to 1/3, and from this expression the bound¬ 
ary conditions for the first GL equation are derived. The 
superconducting current density, defined as the negative 
functional derivative of the GL free energy density with 
respect to the vector potential, for chiral p-wave super¬ 
conductors is 
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= M — 
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where i, j form the canonical base in Gartesian coordi¬ 
nates. The set of operators {a± and az) act on 'ip±, while 
Sy acts on {qj}. The superconducting current contains 
mainly three contributions defined by the following fac¬ 
tors, {K + fci)/4, (^2 + k-i)l2\p2., and (fc 2 — k-i)l^- The 
first one arises from the conventional term in Eq. ([^, 
the second one (we name chiral) is due to the internal de¬ 
gree of freedom (chirality) that appears in Eq. ([^ in the 
form of two nondiagonal terms. Einally, the third contri¬ 
bution arises from the diagonal terms Dy\) in Eq. 

([^, and resembles the Zeeman interaction between the 
magnetic field and the electron spin. In this case ^2 — k^ 
plays the role of the magnetic moment. The boundary 
conditions imposed on Eq. ([^ for our square mesoscopic 
sample are given as: 

at north and south sides, 

at east and west sides. (8) 

It is straightforward to show that the boundary condi¬ 
tions of Eq. ([^ set the perpendicular current at the 
edges to zero, i.e. they impo se spec ular reflection in the 

chiral p-wave superconductorjg important to re¬ 
mark also that they are parameter-independent, so they 
provide the proper boundary conditions for Eq. ([^ but 
also for the most general case of Eq. ([^. With Eq. Q 
we have completed the set of equations needed for the 
GL description of a chiral p-wave mesoscopic supercon¬ 
ductor. Eq. ([^ is numerically solved using finite differ¬ 
ences and the link variables technique of Ref. [15] on a 
square lattice with mesh grid = hy = 0.1. On the 
other hand, the temporal derivative is discretized using 
the Runge-Kutta method of first order. Before conclud¬ 
ing this section, we give the reduced expression for the 
dimensionless free energy, since it allows us to find not 
only the lowest energy (ground) states but also the stable 


^/’+ 0 1 

Dyll}+ + Dy'ijj_ = 0 j 

^/’+ + V'- = 0 1 

0^1^+ - = 0 j 


states with slightly higher energies (metastable states) 
The free energy reads: 

+ + r(^*d,^)2}, (9) 

where Fq = is the bulk free energy at zero field. 


III. ISOTROPIC CASE (CYLINDRICAL FERMI 
SURFACE) 

The results obtained using Eq. Q for a square 8^ x 8^ 
sample are summarized in Eig. showing the dimension¬ 
less free energy and the vorticity of vector 4/ as a function 
of the external magnetic field H. Panels (b) and (c) show 
the vorticity of the ground states of our superconducting 
sample, labeled a-j in Eig. [^(a), where v+{-) is the vor¬ 
ticity of component Note that both and v_ 

remain constant along the stability curves of each state 
in panel (a), and as such are good identification numbers 
for these states. Contour plots in Fig. [^show the order 
parameter 'k corresponding to the ground states a - d. 
While the left and central columns of Fig. [^show contour 
plots of the superconducting density of each component, 
IV'-i-P and IV'-P, respectively, the third column shows the 
difference between the angular phases of the components, 
i.e. 6+ — 9-. 

The ground state a of Fig. shows one anisotropic 
vortex in each component, i.e. vorticity = —1 in com¬ 
ponent ■)/+ and V- = 1 in component ip-. The contour 
plots of the ground state b in Fig. show the vortex free 
state in component ip^ and the giant vortejP^ with vor¬ 
ticity = 2 in component ip-. The subsequent ground 
state c has vorticity i/+ = 2 and i/- = 4, where |V’-P 
contains four vortices close to the corners, meanwhile 
{ipY shows a pronounced depletion around the center 
of the sample. The corresponding phase difference figure 













4 



H/Ho2 



FIG. 1. (Color online) (a) Free energy in units of the bulk 
condensation energy at zero field (Fo) as a function of the 
external magnetic field in units of the bulk upper critical field 
{Hc 2 ), for a square mesoscopic sample of size 8^x8^. Letter 
labels denote different found ground states. Some metastable 
states (not labeled) are also shown in this figure. Vorticity of 
components ij}+ s^nd ip- of the ground states of panel (a) are 
shown in (b) and (c) respectively. The difference in vorticity 
(y+—V- =2) between the components is in perfect agreement 
with the analytically predicted solution ^ ■ 


reveals that the depletion in component is a conse¬ 
quence of two vortices and two vortex-antivortex pairs 
there. The ground state d has six vortices in |'0_P in 
full agreement with the vorticity reported in Figs. ^ (b) 
and (c) (i^+ = 4 and V- = 6 ). However, the density 
fails to convincingly show any signature of a vortex. The 
vorticity zz_|_ = 4 of component '0+ is visible in the phase 
difference figure 2 (d), where 10 discontinuities are found 
along the edges as a consequence of six vortices from ip_ 
and four from Four vortex-antivortex pairs at the 
center of the sample are also visible in this contour plot, 
but do not affect the total vorticity. 

From the comparison between Figs. (c) and (d) one 
sees that with increasing the magnetic field the compo¬ 
nent ip- dominates its partner component '0+. The dom¬ 
inance of Ip- over ' 0 _|_, especially at high fields impedes 
the proper description of the vortex configuration in the 
latter component. In order to describe the components 


FIG. 2. (Golor online) Ground states a - d of Fig. Left and 
central columns show the contour plots of the superconduct¬ 
ing densities components \ip+\'^ and IV'-P, respectively. Right 
column shows the difference between the angular phases of the 
components, i.e. 6+ — 9-. 


of the order parameter on an equal footing, a more suit¬ 
able representation is in terms of ipx and ipy. Figs. [Sland 
l^show contour plots of |' 0 xP, IV'wP and cos( 6 * 3 , — ^ for 
ground states a - j of Fig. Fig.|^(a) shows cos{9x — dy) 
for ground state a (from now on called the phase differ¬ 
ence figure), and reveals a linear domain wall. Its exten¬ 
sion across the sample coincides with the stripe where 
density \'4>y\'^ vanishes. On the other hand, the partner 
component, tpx, is free of vortices. Ground states b and c 
look similar in both densities, although from the compar¬ 
ison between their phase difference figures in Fig. i(b) 
and (c) respectively, we see four domain walls in state c 
and none in state b. The domain walls (DWs) of ground 
state c define a path where the difference between the 
angular phases of components ipx and ijjy are 0 or tt, i.e. 
dx — 9y = 0, TT. Ground state d shows two vortices in 
density |' 0 a;P and none in IV'yPj while its corresponding 
phase difference figure shows the four domain walls of 
ground state c plus two other alternating domain walls 
that weakly connect the former ones. The contour plots 
of Fig. (e), for state e show clearly two vortices in 
each component. They look indistinguishable just from 
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the analysis of their densities, but their phase difference 
figure reveals that there are two vortices, one in each 
component, that combine to produce a different signa¬ 
ture from the remaining vortices. While the uncorrelated 
vortices lead to the formation of the alternating domain 
walls towards sample edges, the pair of correlated vor¬ 
tices align their cores and do not show any domain wall 
between them. The alternating domain wall is therefore 
the signature of the half quantum vortex (HQV), in con¬ 
trast to the other signature without domain wall that 
corresponds to the full vortex (FV). 



FIG. 3. (Color online) Ground states a - e of Fig. plotted 
correspondingly in panels (a) - (e). Left and central columns 
show the contour plots of the superconducting densities com¬ 
ponents I'i/’aiP and I'i/’i/Pj respectively. Right column shows 
cos{6x — &v), where 6x,y are the angular phases of compo¬ 
nents tpx and ipy. 

Fig. a shows the remaining ground states f - j of Fig. 
a Both densities in ground state / clearly show two vor- 



FIG. 4. (Color online) Ground states f - j of Fig. ^ plotted 
correspondingly in panels (f) - (j). Displayed quantities are 
the same as in Fig. 


tices in each component, which are indeed four HQ Vs 
according to the corresponding phase difference figure. 
Ground states g and h show one common feature, hav¬ 
ing different number of HQVs per component, but all of 
them aligned vertically in component il’x and horizontally 
in component ij^y. On the other hand, the corresponding 
phase difference figures for states g and h show that: (i) 
two vortices, one per component, combine to form one 
FV in state g, and (ii) four vortices, two per each com¬ 
ponent, combine to form one skyrmion in state h. The 
signature of the skyrmion is shown here for the first time: 
four alternating d omai n walls which are connected into a 
circular structure! ^° * ^^ * The skyrmion state here of course 
differs from those of magneti c mater ials due in physics 
and the formation mechanism. F^ I ^^I^ ^ Nevertheless, their 
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topological properties remain similar, as will be presented 
later. The phase difference figure of the ground state i 
shows four DWs around the corners, four HQ Vs close 
to the edges and three FVs in the center. What draws 
attention in all three contour plots of Fig. (i) is that 
there are five vortices in each component (fractional vor¬ 
tices), and among them three align their cores to form 
FVs according to the corresponding phase difference fig¬ 
ure. The triangular array formed by them resembles the 
consequences of vortex-vortex repulsion in conventional 
type II superconductors. Therefore, this supports our 
initial premise that the FV in our analysis is the usual 
Abrikosov vortex of conventional superconductivity. Fi¬ 
nally, the phase difference figure of the ground state j 
shows four DWs, six HQVs and two FVs. One system¬ 
atic comparison of the phase difference figure of ground 
states / - j clearly shows that HQV and FV are indeed 
very different states. While FVs are formed in the sam¬ 
ple center, being favored by confinement, all the HQVs 
remain close to the sample edges. In order to explain 
this difference the following subsection discusses the cal¬ 
culated superconducting currents in the sample. 



a:/? 


FIG. 5. (Color online) Superconducting currents around: (a) 
the upper right DW of Fig. I^c), (b) the upper HQV of Fig. 
3 'd), (c) the FV of Fig. and (d) the skyrmion of Fig. 

4». 


So far, DWs, HQVs, FVs and skyrmions have been 
distinguished in this work according to their signatures 
in the phase difference plots. The superconducting cur¬ 
rent, the physical quantity intertwined with the magnetic 
field, also allows us to identify more characteristic fea¬ 
tures of the novel topological solutions. Figs. i (a)-(d) 
show the supercurrents around one DW, HQV, FV and 
skyrmion, respectively. Fig. [^(a) zooms in the supercur¬ 
rents around the right-top DW of Fig. [^c). Its super- 



FIG. 6. (Color online) Diagonal profiles of the contour plots 
IV'iP corresponding to ground states c and h of Fig. shown 
in panels (a) and (b), respectively. Blue and green arrows 
indicate the DW and vortex core locations, respectively. 
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FIG. 7. (Color online) Contour plots of the magnetic induc¬ 
tion corresponding to: (a) the supercurrents of the DW of 
Fig. [^a), (b) the HQV of Fig. [5];b), the FV of Fig. |^c), and 
(d) the skyrmion of Fig. id). 


currents are flowing from the upper left to the bottom 
right corner, where the DW outer and inner screening 
currents flow counter-clockwise. Aside the DW current, 
there are two narrow stripes where the supercurrent van¬ 
ishes, consequence of the continuity of the supercurrent 
field and the compensation between the screening and the 
DW currents. In order to understand better the origin 
of the DW currents. Fig. [^a) shows the line profiles of 
the corresponding superconducting densities |V’±P along 
the diagonal line defined by y = a:. Light (green) ar¬ 
rows point towards the already seen vortex cores of com¬ 
ponent ip- in Fig. [^c). Dark (blue) arrows indicate 
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FIG. 8. (Color online) Textures of the ground states h and g of Fig. according to the mapping n = where a 

are the Pauli matrices. Colors show the amplitude of the z-component of h. 


the center of two DWs defined by the intersection where 
the densities IV'-P and |'0+P become equal. A straight¬ 
forward calculation that exploits the previous definition: 
IV'-i-P = IV'-P) leads to: — V'yV'S = 0, or equiva¬ 

lently: sin (Ox — Oy) = 0, in full agreement with the con¬ 
tour plots of the phase difference in Fig. [^c), where 
dx — 0y = O, or 9x — 0y = 'n along the DWs. 

The magnetic induction that corresponds to the DW 
supercurrents of Fig. j^a) is shown in panel (a) of Fig. 
[Tj It is calculated using the Maxwell equation, 

k‘^VxB = J, ( 10 ) 

where = K^jd, with k = 46 being the GL parame¬ 
ter reported for SRO along the c axis,El and d being the 
sample thickness which we suitably choose to be 2^. The 
contour plot of Fig. [^a) shows that the magnetic in¬ 
duction corresponding to the DW is weak and strongly 
screened by the Meissner effect. This fact represents an 
obstacle for the detection of DWs signatures in direct 
measurements of their magnetic response such as in mag¬ 
netic force microscopy (MFM) or scanning Hall probe 
microscopy (SHPM). 

Fig. I^b) zooms in the supercurrents around the up¬ 
per HQV of Fig. [^d). It shows two adjacent counter¬ 
flowing streams with the bottom one flowing clockwise 
and belonging to the HQV supercurrents, while the top 
one flows counter-clockwise and represents the screening 
currents. The Meissner effect for the HQV is anisotropic 


due to the boundary conditions of Eq. (§. From the 
supercurrent equation ([^, and the local approximation 
tpy « 0, or drawn from Fig. |^d), one eas¬ 

ily obtains: J « Im{ipP^Dx4’+i + ^'4’+Dy'4’+j}- After 
straightforward calculations and replacing the covariant 
derivative once again one obtains: J « \tp+\‘^[{dxd i + 
\dy9j) + — ^coscpj)], which draws attention 

since the screening currents are defining elliptical equipo- 
tential lines. Thus, the anisotropic screening of the super¬ 
conductor towards the HQVs causes them to move along 
the easy-screening direction which in this case is along y. 
The contour plot of the magnetic induction correspond¬ 
ing to the supercurrents of Fig. [^b) is shown in Fig. 
Qb). As expected from the two counter-flowing streams 
seen in the HQV supercurrents, the magnetic induction 
also shows adjacent local maximum and local minimum. 

Fig. i c) zooms in the superconducting currents 
around the FV of Fig. (g), and shows that the FV cur¬ 
rents flow clockwise and vanish as we move away from the 
FV core. This vanishing is due to the spatially isotropic 
Meissner effect, unlike in a HQV, that screens the FV 
currents. As expected, its magnetic induction signature 
[see Fig. [^c)] agrees well with that of the Abrikosov 
vortex. 

The superconducting currents around the skyrmion of 
Fig. i (h) are shown in Fig. (d). Unlike the FV, 
the skyrmion supercurrents clearly show outer and in- 







































ner structures. The supercurrents of the outer struc¬ 
ture flow clockwise while the supercurrents of the inner 
flow counter-clockwise. The skyrmionic DW of Fig. Sh) 
along with its supercurrents in Fig. id) shows cylindri¬ 
cal symmetry, and one easily deduces that the same sym¬ 
metry is present in densities |V’±P- Line profiles of |V’±P 
then provide enough information to unveil the skyrmion 
supercurrents [see Fig. (b)]. The inner structure of the 
skyrmion is defined by: '0- =0 and V0+ =0, where 0+ is 
the angular phase of component 0+, which when replaced 
in the supercurrent of Eq. Q yields J « 
with 0 being the azimuthal vector in cylindrical coordi¬ 
nates. Hence, the supercurrents at the skyrmion core are 
mainly due to the Meissner effect. However, away from 
the skyrmion core the scenario changes since the circular 
DW of Fig. 1^ (h) is met, as indicated by arrows in Fig. 
(b). Close beyond the circular DW, we find that compo¬ 
nent 0+ drops to zero. Replacing in Eq. 0 0+ = 0 and 
bearing in mind that one giant vortex is hosted in com¬ 
ponent 0_, the supercurrent in cylindrical coordinates 
becomes: J « |0_p(— ^ -I- ^r))(j). The magnetic 

induction corresponding to the skyrmionic supercurrents 
of Fig. [^d) is shown in the contour plot of Fig. [^d). 
It clearly shows one local minimum at the skyrmion core 
surrounded by one circular stripe of local maxima, and as 
such can be directly imaged in magnetic measurements. 

Finally, before concluding this section, we briefly de¬ 
scribe the superconducting order paramet er in terms 
of the 3-dim real vector field n, defined byfi^Illl 


n = 




( 11 ) 


which maps the complex spaces C x C of components 
0a;, 0J, into the real space K^. A straightforward cal¬ 
culation yields: h = sinQ;cos0, sinasin^, cosa, where 
sin a — cQg Q, — -Ibal „ j (h — 0 — 0 

Then, the target space of mapping 0 is the 2-dim 
sphere of radius one, topological invariant 

of the spac es tha t result from mapping ( [Tl] ) is defined by 
the integralUH^l 


Q = ^ [h ■ {dxfi X dyfi) dx dy, (12) 

47r J 

which is widely known as the Hopf invariant. One conve¬ 
nient interpretation of this topological invariant is that it 
counts the number of times that the 3-dim real held {n) 
wraps around the 2-dim sphere (S^). Left and right pan¬ 
els of Fig. [^show the texture fi for ground states h and g 
of Fig. 1^ respectively. The texture for the skyrmion (left 
panel) differs from the texture for the FV (right panel) 
owing to the alternating circular DW characteristic of the 
former state. While at the skyrmion core, field n points 
towards —j, outside the skyrmion it points towards j. 
Along the DW that separates the skyrmion core from 
the outside, the field texture whirls, therefore providing 
to the space the topological charge Q = —2. The field 


texture that corresponds to the FV shows four lobes C4 
symmetric profile where n changes smoothly. Unlike the 
skyrmion and in agreement with our earlier results, the 
field texture for the FV does not show any signature of 
a domain wall separating unequivalent outer and inner 
regions. 


IV. FIELD-DRIVEN TRANSITIONS BETWEEN 
SKYRMIONIC AND VORTICAL STATES 


In bulk and type H superconducting samples vortices 
with phase windings higher than 2n7r, where n is integer, 
are energetically disfavored. The superconductor prefers 
two distant vortices each with phase winding 27r rather 
than one single vortex (giant vortex) with phase winding 
Att. Nevertheless, in samples with dimensions of the or¬ 
der of the superconducting coherence length (mesoscopic 
samples), giant vortices can appear as stable configura¬ 
tions. The stabilization is provided mainly by the con¬ 
finement due to the small sample size, although the exter¬ 
nal magnetic field also contributes through the screening 
currents and the confining force they exert on vortices. 
Field driven transitions from states with multiple di stant 
vortices to giant vortices have been widely reported !^^ * ^^ ! 

In this work we first report the field-driven transitions 
from HQV to FV states. Fig. (a) shows the energy of 
state / of Fig. (a), along with some of its neighbor¬ 
ing states. Panel (b) shows the second order derivative 
of the energy with respect to the external field only for 
state / . While the energy of state / is continuous, its 
second order derivative shows discontinuities indicating 
transitions between distinct states. Three different states 
can be easily distinguished, which we labeled by a circle, 
square and triangle marker. The corresponding distribu¬ 
tions of the superconducting order parameters are also 
shown in the figure: logarithmic contour plots of |0a;P 
and |0j,P are shown in the left and central columns, while 
the cosine of the phase difference is shown in the right 
column. State (Q) shows two fractional vortices in each 
component rendering four HQ Vs according to the phase 
difference contour plot. State (□) shows two HQVs and 
two FVs. The FVs are composed of two fractional vor¬ 
tices belonging separately to each component. The frac¬ 
tional vortices composing the FVs are slightly misaligned 
as can be seen in the density figures. This makes the FVs 
display a small closed domain wall in the phase difference 
contour plot. At high fields the screening currents con¬ 
fine even more the superconducting configuration of state 
(□) transforming it into one state with three HQVs and 
one FV (A). Due to the strong screening currents the 
upper FV of state (□) loses one of its fractional vortices 
which renders one HQV in state (A). The strong confine¬ 
ment also forces the alignment of the fractional vortices 
composing the FV of state (A). 

Another field-driven transition from skyrmion to FV 
state is presented in Fig. 10 Panel (a) shows the energy 
of state j of Fig. along with some of its neighboring 
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FIG. 9. (Color online) Field-driven transition from HQV to 
FV due to confinement in a square mesoscopic sample of size 
8^x8^. (a) Energy of the state / ofFig. [^(a), along with some 
of its neighboring states, (b) Second order derivative of the 
energy with respect to the external field showing three distinct 
states indicated by circular, squared and triangular symbols. 
The corresponding components of the superconducting order 
parameter are shown in panels (0)i (D) and (A). Displayed 
quantities are logarithmic contour plots of \'>px\^ and \ipy\^ in 
left and central columns, respectively, while the cosine of the 
phase difference is shown at the right column. 


states. Panel (b) shows the first and second order deriva¬ 
tives of the energy with respect to the external field only 
for state j. Unlike in Fig. [^(b), here the second derivative 
is continuous as well as the first derivative. Nevertheless, 
this does not mean that there are no distinct states along 
the stability curve of state j. Circle, square and triangle 
markers (Q: ^ and A) indicate three states at weak, 
intermediate and strong confinement, respectively. At 
weak confinement the phase difference figure shows six 
HQVs and one skyrmion (see Fig. 10 (O))- interme¬ 
diate confinement, state (□) shows in |'0a;P that two out 
of the four fractional vortices composing the skyrmion 
of panel (Q) have merged into one single discontinuity. 


FIG. 10. (Golor online) Another example of a field-driven 
transition between skyrmionic and vortical states along the 
state j of Fig. Displayed quantities are the same as in Fig. 
1^ with the only exception that in panel (b) the first order 
derivative of the energy with respect to the external field is 
also shown. 


This merger of initially distant fractional vortices renders 
the domain wall of the skyrmion asymmetric. At strong 
confinement (A) the former fractional vortices split their 
cores along the horizontal axis. According to the phase 
difference figure they join two other fractional vortices in 
density IV'yP to form two horizontal FVs in the center of 
the sample. As can be easily seen, the vorticity of the 
superconducting components along this field-driven tran¬ 
sition is constant, unlike in Fig. where it was not. This 
fact explains why the second order derivative is continu¬ 
ous here and discontinuous in Fig. 


V. TEMPORAL DYNAMIC TRANSITIONS 

To date, no works have treated the time-dependent 
phenomena within the GL formalism for chiral p-wave 
superconductors. Here we benefit from the temporal 






















































10 


evolution included in the TDGL equations to report for 
the first time dynamic transitions involving vortices and 
skyrmions. 



FIG. 11. (Color online) Temporal vortex-skyrmion transition 
in a square mesoscopic sample of size 12^x12^. Panel (a) 
shows the free energy of states i and j containing 10 and 12 
fractional vortices per component, respectively. The energy 
of state i is discontinuous a.t H ~ 1.06iFc2 reflecting a first 
order transition. Panel (b) shows the temporal evolution of 
the energy at the latter transition. Three states, initial, inter¬ 
mediate and final are denoted by circle, square and triangle 
markers, respectively. The components of the superconduct¬ 
ing order parameter corresponding to each state are shown in 
panels (0)i (D) and (A). 

The dimensionless free energy as a function of the ex¬ 
ternal field for states i and j is shown in panel (a) of 
Fig.[TTl Unlike in Fig. [^the sample size here is 12 ^x 12 ^ 
rather than the 8 ^x 8 ^, which was a suitable choice to 
study the evolution of the superconducting configuration. 
Panel (b) shows the temporal evolution of the free energy 
at the discontinuous step in energy in panel (a). Three 
states, initial, intermediate and final, are denoted by cir¬ 
cle, square and triangle markers. The corresponding su¬ 
perconducting order parameters are shown in panels (0 )j 


(□) and (A), respectively. The displayed quantities in 
the latter panels are the same as in those of Fig. The 
initial state (Q) is a multi-vortex-skyrmion state contain¬ 
ing two pairs of skyrmions and FVs, surrounded by eight 
HQVs at the sample edges. This state was not obtained 
for sample size 8 ^ x 8 ^ mainly due to the strong confine¬ 
ment there. At the intermediate state (□) two fractional 
vortices nucleate in each component of the superconduct¬ 
ing order parameter forming two FVs according to the 
phase difference contour plot. The four FVs of the inter¬ 
mediate state then combine following the inverse process 
of the one described in Fig. |10[ to form two skyrmions 
as depicted in state (A). In the Supplementary Mate¬ 
rial, we present the animated data showing the temporal 
evolution of the superconducting order parameter, as a 
more convenient view of the transition reported here. It 
is noteworthy here that all field-driven transition from 
HQV or skyrmion to FV states and vice-versa are essen¬ 
tially driven by HQV penetration and recombination into 
other topological entities. 


VI. ANISOTROPIC CASE 


A. Strong chiral limit 


This far, the ground states of a p-wave mesoscopic su¬ 
perconductor with size 8 ^ x 8 ^ have been obtained under 
the assumption of weak coupling and with a cylindrical 
Fermi surface, which led us to set the ki parameters to 
1/3. However, several works have reported or suggested 
other scenarios for SRO such as: (i) multiband super¬ 
conductivity wit h the ID Fermi sheets developing super¬ 
conducting order, 12 Sl^ or (ii) anisotropy in the cylindri¬ 
cal Fermi surface.^ In order to include just anisotropy 
in the Fermi surface, while preserving single-band super¬ 
conductivity, in this section we introduce the parame¬ 
ter {6k), which sets the ki’s to: ki = 1/3 + 26k, and 
^2 = fca = 1/3 — i5fc. The motivation behind this choice 
is that the theoretical values for the k^ parameters cor¬ 
responding to the three Fermi sheets ( 7 , a and /3) lie 
between 1/3 < fc/ < 1 and 0 < /c“’^ < 1/3,1^^ respec¬ 
tively. The GL equation for p -wave superconductors with 
anisotropy in the Fermi surface becomes: 






+ 6k[& - 2 (n 2 d+ + n^d_)]^. 


(13) 


By tuning 6k within the interval [0,1/3], the strength of 
the chiral terms is changed, therefore driving the system 
between two limiting cases: the left limiting case being 
at = 0 and given by Eq. (§, and the right limiting 
case being aA, 6k = 1/3 where the chiral coupling between 
the superconducting components is set to zero. 

Fig. [ 1 ^ summarizes the results obtained from the 
simulations that numerically approach Eq. (131 with 


6k = 0.03. The energy against field plot of Fig. 12 
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FIG. 12. (Color online) The free energy as a function of the 
external magnetic field, showing ground states b - j plus one 
met astable state a, from the numerical simulations using Eq. 
| |13| l with Sk = 0.03. The parameters kt thus only slightly 
deviate from the value 1/3 obtained when a cylindrical Fermi 
surface is considered. Panels (a) and (b) show the supercon¬ 
ducting density components I'i/’+P I'l/’-P of the states a 

and b, respectively. 


there were obtained a.t H = 0, thus those are the sponta¬ 
neous currents widely sought in experiments. The spon¬ 
taneous currents are composed mainly of two counter¬ 
flowing streams at left and right sides of our sample. 
They are the chiral edge currents predicted by Furusaki 
et. aZ.j^Matsumoto and Sigrist,l®and Stone and Roy.^^ 
Along the line x = 4^, the linear domain wall (DW) of 
Fig. [^a) separates the left and right sides showing an en¬ 
hancement in the supercurrents around the center. The 
magnetic induction corresponding to the supercurrents 
of panel (a) is shown in panel (b) of the same figure. 

Panel (c) of Fig. 13 shows line profiles of the magnetic 
induction and the y-component of J along the line y = 
4^. This plot agrees well with the result of Matsumoto 
and Sigrist which showed that Jy {Bz) is an even (odd 
function of x along the line perpendicular to the DW.b 
Finally, panels (d) and (e) provide important information 
that allow us to calculate the supercurrent along the DW. 
From panel (d) the DW is defined by \'4>+\ = IV"-! at 
X = 4^, but along this line panel (e) tells us that not only 
the magnitudes of the superconducting components are 
equal but also their angular phases. Then, from Eq. 0 
our estimation for the superconducting current along the 
linear domain wall is simply Jy{x = 4^) = ki\ip+\’^dy9+. 


shows nine ground states labeled by letters. The com¬ 
parison between Figs. and reveals one important 
fact: the energy of the state a is higher than the energy 
of its adjacent state b. Actually, state a here is no longer 
the ground state at low fields i/ « 0, unlike in Fig. 
where it was. Contour plots of the superconducting or¬ 


der parameter (4/) that correspond to the states a and b 
of Fig. 12 are depicted in insets (a) and (b), respectively. 
The comparison between the insets of Fig. and the 
corresponding states in Fig. shows that despite of the 
small anisotropy introduced in the GL equation, the su¬ 
perconducting configuration of these states is practically 
identical in both cases. 


Two decades have passed since the discovery of the 
unconventional properties of strontium ruthenate, but to 
date there has not been a conse nsus w hether or not it 
is a chiral p-wave superconductorThe main experi¬ 
mental results that support unconventional superconduc¬ 
tivity in SRO are provided by the set of measurements 
carried out using t echni ques such as the nuclear magnetic 
resonance (NMR) 12213^ /rSR,l^ the Kerr effect,and can¬ 
tilever magnetometry.l2SI The smoking gun evidence that 
lacks, and which, if found, would convince the scientific 
community is the finding of the t heoretic ally predicted 
spontaneous currents in Sr 2 Ru 04 . h^F 6 | 37 | interestingly, 
what we just found in this work is that the state with 
spontaneous currents is no longer the ground state when 
the GL model slightly deviates from the isotropic case at 
iL « 0, i.e slightly deviated from the cylindrical Fermi 
surface. This energy lift of the state with spontaneous 
currents makes it even harder to be detected. Fig. [T^a) 
shows the supercurrent distribution corresponding to the 
state a of Fig. We note that the currents displayed 


B. Strong Zeeman limit 


Microscopy with superconducting quantum interfer¬ 
ence devices (SQUIDs) and scanning Hall probes (SHPs) 
have recently detected vort ex coa lescence in single crys¬ 
tals of strontium ruthenate.l^SEH One possible explana¬ 
tion for this behavior is the existence of at least two dif¬ 
ferent coherence lengths arising from multigap supercon¬ 
ductivity, and which lead to attractive (repulsive) inter¬ 
action at long (short) ranges.l^ Refs. and [31] have 
reported that within their corresponding resolutions no 
convincing evidence for spontaneous currents and DWs 
has been found yet. In order to explore more super¬ 
conducting configurations, comprising DWs, HQVs, FVs 
and skyrmions as the fundamental entities, in what fol¬ 
lows a different set of the ki parameters is defined by: 
ki = 1/3, ^2 = 1/3-1-i5fc and kz = 1/3 —Such a choice 
of parameters enables one to keep constant the strength 
of the chiral terms while varying 5k. The first GL equa¬ 
tion for this particular choice of parameters reads 


dt^ = \[ d ^ + n^<7+ + n^d_ - ^Haz 


4 - 


+ ^(l- 


3 | 4/|2 ^ 

4 4 /■ 


(14) 


The resemblance of the fourth term in the upper line 
of Eq. (141 with the Zeeman term leads us to interpret 
5k here as a sort of magnetic moment. In order to study 
the dependence of the superconducting configuration on 
the anisotropy parameter, in Eq. ( |14[) th e magnetic field 
is kept fixed while 5k is varied. Fi g. |14[ a) plots the free 
energy of the states, solving Eq. (|14[), as a function of 
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FIG. 13. (Color online) (a) Superconducting currents corresponding to the state a of Fig. 12 These currents, which were 


obtained at zero field, are composed of two edge currents with different chiralities and flowing in opposite senses, (b) Contour 
plot of the magnetic induction {Bz) calculated from the supercurrent s of panel (a), (c) Line profiles of Jy and Bz along the line 


y = 4^. (d) Line profiles of \'4’±\'^ corresponding to the state a of Fig. 
'il)± along the line a: = 4^. 
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(e) Line proHles of the angular phases of components 


the anisotropy parameter 6k- Circle, square and triangle 
markers denote three states whose IV'+P and IV'-P diag¬ 
onal {y = x) line profiles are shown in panels (b) and (c), 
respectively. From panel (c), and unlike in panel (b), one 
clearly sees that for high values of 5k the density |V’-P 
diminishes. Our explanation for this behavior is through 
the definition of two effective coherence lengths, one for 
each superconducting component. Defining them as the 
coefficients in front of the linear terms V’+ and '0_ in 
Eq. (141, they read: = 1 — 6kH and ^_ = 1 -I- 5kH, 

respectively. With H fixed and 5k increasing, (f_) 
becomes smaller (larger) therefore leading to an effec¬ 
tive reduction (increase) of confinement in component '0+ 
(V'-). Concerning the phase, contour plots of 0+, 9- and 
cos {9x — 9y) corresponding to the states denoted by cir¬ 
cle, square and triangle markers are shown in rows (0)j 
(□), and (A). According to the phase difference figure, 
state (O) is composed of two concentric skyrmions, one 
circular and one rhomboidal. From the contour plot of 
9_ one sees that the circular skyrmion arises from the 
formation of one giant vortex in 0_ with phase wind¬ 
ing drr. The phase difference figure corresponding to 
state (□) shows one irregular closed domain wall emerg¬ 
ing from the intersection of the circular and rhomboidal 
skyrmions. Its formation is determined by the annihi¬ 
lation of the giant vortex in that has split into two 
fractional vortices. Finally, the phase difference figure of 
state (A) shows four FVs with cores slightly asymmet¬ 
ric as can be seen from the small circular DWs present 


there. Due to the density |0-P has been substantially 
depleted at this value of dk, the superconducting state is 
completely defined by component il)+. Hence, what we 
have achieved by adding anisotropy to the chiral p -wave 
model of Eq. (§ , is a chiral polarization enhanced due to 
the strong confinement present in a mesoscopic sample. 


VII. CONCLUSIONS 


In summary, we have studied in detail the 
Ginzburg-Landau m odel that describes p -wave tetrag¬ 
onal superconductor^^E^ and all the possible states of 
a mesoscopic superconducting sample as a function of 
the external magnetic field and the anisotropy param¬ 
eters of the material. Due to the odd parity the or¬ 
der parameter is a two-component complex vector and 
the fundamental solutions of the corresponding TDGL 
equations, that we obtained numerically, are fractional 
vortices, i.e. solutions where the phase winding 27r is 
found in one component but not in the other one. In 
two- and three-band superconductors similar fractional 
vortices were obtained between components, but for dif¬ 
ferent reasons .EsHUl Fractional vortices in different com¬ 
ponents can combine to form a cored/full-vortex state, 
as well as a coreless/skyrmion state seen in phase dif¬ 
ference and magnetic response figures. Skyrmions arise 
when same number of fractional vortices in each compo¬ 
nent combine to form a closed domain wall that sepa- 
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FIG. 14. (Color online) (a) Free energy as a function of the anisotropy parameter 5k with the external magnetic field fixed at 
H = 0.530[iyc2]. Three distinct states are indicated by circular, square and triangular markers (Oi D and A), (b) and (c) show 
line profiles of I'i/’+P a^nd respectively, along the diagonal line {y = x) corresponding to the states of panel (a). Columns 

(0)i (n) and (A) show contour plots of and cos — Qy) corresponding to the denoted states of panel (a). 


rates d istinc t intercomponent phase difference {6x — Oy) 
regions .1221211 Alternating segments between 1 and —1 in 
the cos {Ox — Oy) between fractional vortices along the do¬ 
main wall is the main signature for skyrmions. While 
for skyrmions the topol ogical charge (Q) is defined by 
the hopf invariant for vortices it is by the circu¬ 
lation of the superconducting velocity. The screening 
currents of half-quantum vortices are anisotropic and in 
Cartesian coordinates we have analytically shown that 
the equipotential lines of the screening currents are el¬ 
lipsoidal rather than circular as in full vortices. This 
anisotropic screening leads half quantum vortex to be 
attracted towards the sample edges. The mesoscopic size 
of our samples provides stability to the half quantum 
vortices and the Q = 2 skyrmions, in contrast to larger 
systems where larger values of Q were considered,^ and 
bulk systems where the half quantum vortic es have been 
usually regarded as high-energy states.l^I^I! Actually the 
mesoscopic size of the sample plays a remarkable role in 
the stability of skyrmions as well as in the here reported 
novel transitions (e.g. from a skyrmion to a full vortex). 
At high external fields, above the Hd critical field, states 
with different configurations of skyrmions and half quan¬ 


tum vortices gradually transform into full vortices owing 
to the increased screening currents and confinement ef¬ 
fects. 

To date, the only superconductor expected to be 
p-wave is strontium ruthenate, with eno ugh e vidence 
demonstrating its unconventional behavior.EEHMl Nev¬ 
ertheless, many works have failed to convincingly de¬ 
tect spontaneous currents, half-quantum vortices and 
skyrmions in large samples. What we have demonstrated 
here is that: (i) even by slight anisotropy in the Fermi 
surface, the state with spontaneous currents is no longer 
the ground state at « 0, and (ii) field-driven tran¬ 
sitions between half-quantum vortex, full vortex, and 
skyrmions provide a better method to indirectly prove 
their existence in magnetic measurements. 
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